Poster Presentation 


International Conference on Clouds and Precipitation 
July 7-1 1, 2008, Cancun, Mexico 


ROLE OF GRAVITY WAVES IN DETERMINING CIRRUS CLOUD 

PROPERTIES 

David O’C. Starr 1 , Tamara Singleton 2 and Ruei-Fong Lin 3,1 

Laboratory for Atmospheres 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 

2 Department of Mathematics 
University of Maryland, College Park 

3 Goddard Earth Sciences and Technology Center 
University of Maryland, Baltimore County 


Abstract 

Cirrus clouds are important in the Earth’s radiation budget. They typically exhibit variable 
physical properties within a given cloud system and from system to system. Ambient vertical 
motion is a key factor in determining the cloud properties in most cases. The obvious exception 
is convectively generated cirrus (anvils), but even in this case, the subsequent cloud evolution is 
strongly influenced by the ambient vertical motion field. It is well know that gravity waves are 
ubiquitous in the atmosphere and occur over a wide range of scales and amplitudes. Moreover, 
researchers have found that inclusion of statistical account of gravity wave effects can markedly 
improve the realism of simulations of persisting large-scale cirrus cloud features. Here, we use a 
1 -dimensional (z) cirrus cloud model, to systematically examine the effects of gravity waves on 
cirrus cloud properties. The model includes a detailed representation of cloud microphysical 
processes (bin microphysics and aerosols) and is run at relatively fine vertical resolution so as to 
adequately resolve nucleation events, and over an extended time span so as to incorporate the 
passage of multiple gravity waves. The prescribed gravity waves “propagate” at 15 ms' 1 , with 
wavelengths from 5 to 100 km, amplitudes range up to 1 m s' 1 . Despite the fact that the net 
gravity wave vertical motion forcing is zero, it will be shown that the bulk cloud properties, e.g., 
vertically-integrated ice water path, can differ quite significantly from simulations without 
gravity waves and that the effects do depend on the wave characteristics. We conclude that 
account of gravity wave effects is important if large-scale models are to generate realistic cirrus 
cloud property climatology (statistics). 
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1. INTRODUCTION 

The deposition rate, % , (in units of 
mass of condensate per mass of air per unit 
time) is an essential component in cirrus 
evolution. If the deposition rate is known, we 
could explain the evolution of the observed 
moisture field in a more quantitative fashion. 
The goal of this research is to present a 
new method for estimating the deposition 
rate using ground-based remote sensing 
measurements. 

2. MEASUREMENTS 

We use measurements of the extinction 
coefficient, b exl L , from a Raman lidar (RL) 
and the equivalent radar reflectivity factor, 
Z e , from a 35 GHz Millimeter-wave Cloud 
Radar (MMCR) to derive x ■ Both 
instruments are located at the ARM Climate 
Research Facility (ACRF) Southern Great 
Plains (SGP) site. The different but 
complimentary operating wavelengths 
(Donovan and van Lammeren, 2001) make 
the retrieval of x possible. Details for both 
the ARM RL and MMCR can be found at 
http://www.arm.gov . 

The parameter b alL is computed from 

the ARM RL nitrogen channel using a 
procedure described in Ansmann et al. 
(1992). The water vapor mixing ratio q v is 

computed using the ratio of the lidar signals 
at 408-nm (water vapor scattering) and at 
387-nm (nitrogen scattering). Temperature 
profiles are interpolated between four times 
a day radiosonde profiles. 


3. METHOD 


In this study, we consider columns and 
bullet rosettes because we are mostly 
interested in the upper tropospheric cirrus 
formed in situ. 

The modified Gamma function has 
been widely used to describe the size 
distribution of ice particles in cirrus clouds, 


i.e., N(D)dD = 


exp( — — ) 


D_.dD_ 
A A 


where D is the particle maximum 
dimension, N, is the number of particles 
per unit volume, v is a parameter 
controlling the width of the distribution and 
A is a characteristic length. Note that 


A(v-l) is the modal length. An 
advantage of this practice is that a simple 
formula results for the moments. Thus, if a 
particle property can be formulated or 
approximated by a power-law function (i.e. 
x = aD p ), its integrated property over the 
entire size distribution can be obtained 
easily. 

The extinction coefficient b^ lL is a 
moment of the particle size distribution 
(PSD), b ext L ^2\A{D)N{D)dD, where A is 

the particle cross section and N{D) is the 
size distribution function. The power-law 
relationship between A and D has been 
observed from various studies and is 
summarized in Heymsfield and Miloshevich 
(2003). 



Analogously, Z e is also a moment of the 
PSD, Z e = - \cr R {D)N(D)dD , 

* I I 

where a R is the backscatter cross section, 
X is the wavelength, and | k w | 2 is the 
dielectric constant for water particles. The 
power-law relationship between a R and D 
has been given in Aydin and Walsh (1999). 

The power law relationship is suitable to 
describe the relationship between the ice 
crystal growth and its maximum dimension 
because the approximation error is within 
10% when using 2 or 3 segments to capture 
the full parameter range. In this work, the 
capacitances of bullet rosettes and columns 
are formulated according to Chiruta and 
Wang (2003) and the spheroid 
approximation, respectively. The actual 
value of the deposition coefficient, , is still 

a debatable question. Thus, we use a 
sensitivity test approach to examine its 
effect. It is set to 0.006 (slow growth, e.g., 
Magee et al., 2006), 0.1, or 1 (fast growth) 
in the calculations. 

The flowchart of the procedure is briefly 
illustrated in Figure 1. First, R , the ratio of 
bexi,L ancl Z e times a constant, is used to 

obtain D, , which is then used to acquire 
N, . The parameter v is set to 2 after tests 

confirm a less than 10% sensitivity to the 
parameter by varying its value between 1 
and 3. Note that in the procedure, the 
incomplete gamma function is used to 
account for segments of power-law 
relationship. A similar approach has been 
adopted to retrieve ice water content (IWC) 
and mass median size using the zeroth and 
first moments of the Doppler spectrum from 
a millimeter-wavelength Doppler radar 
(Mace et al., 2002). 

4. RESULTS 

A cirrus observed on 7 December 1999 
is examined in this study. A modeling study 
of the same cloud system is reported in 
Comstock et al. (2008). The RL and MMCR 


measurements are plotted in Figure 2. Note 
that the uncertainty for the RL measured q v 

increases from ~10% at z = 8.5 km to ~30% 
near the cloud top. Roughly 25% of cloudy 
pixels in the upper cirrus are highly ice 
supersaturated (RHI>120%). 

Spotty high extinction (6 w<i ) regions are 

located near the cloud top. From our 
calculations, these bright spots contain 
many small particles (Figure 3). 
Nevertheless, 80% the cloudy pixels contain 
N \ between 5 and 500 L' 1 where only 10% 

of the cloudy pixels contain TV, greater than 
500 L' 1 . 

Also shown in Figure 3 is a decreasing 
trend of D , with respect to height, indicating 
the importance of sedimentation-growth in 
this cloud. 90% the cloudy pixels contains 
PSDs of D, less than 100 micron. 

X in Figure 3 is calculated assuming 
column particles and /?, = 1, which is the 

fastest growth scenario. The results indicate 
that, if those assumptions on the particle 
shape and the deposition coefficient are 
adequate, particles are uptaking water 
vapor very efficiently in those high ice 
supersaturation regions. However, this 
seems to directly contradict the frequent 
(25%) and non-spotty occurrence of high ice 
supersaturation in this cloud. Does this 
contradiction imply that should be 

smaller than 1? It is yet too premature to 
make such a conclusion. More studies are 
needed to understand this issue. 

One naturally wonders whether the 
significance of latent heat release in the 
cirrus evolution might change with /?,. . From 

Figure 3, persistent high positive deposition 
rate occurs in the upper cloud between 
0300 and 0400 UTC. The hourly averaged 
latent heating rate profiles computed from 
six scenarios for this period of time are 
compared in Figure 4. As the figure 
indicates, the differences are enormous. 
Latent heat release in the fast growth 
scenario would contribute significantly to 
buoyancy production and local circulation 



whereas it plays a minor role in the cloud 
dynamics in the slow growth scenario. 

The last panel in Figure 3 shows that a 
moderate (as compared to the local 
moisture tendency) magnitude of 
sublimation takes place near the cloud 
base. Does particle sedimentation and local 
sublimation result in the downward 
development of the cloud observed by the 
ground-based instruments? To answer this 
question, we need to examine the local q v 

tendency first. Most of the pixels below z = 
8.2 km has a random error less than 10%. 
According to our estimate, the local 
tendency calculated from two consecutive 
pixels in time has an error less than 20% 
when their difference is greater than 20%. 
We’ve found the uncertainty in the 
measured q v does not hamper our analysis 

because this cloud base experiences a fast 
moistening that is within the measurement 
limit. According to our calculations, even the 
fastest growth scenario is not able to 
explain the drastic moistening near the 
cloud base. The moistening should have 
been a result of sub-grid (relative to the 
large-scale data derived from the 
constrained variational analysis) q v forcing, 

i.e., likely the unresolved horizontal 
advection. 

5. SUMMARY 

A new method that estimates the 
deposition rate using data collected from 
Raman lidar and the millimeter wave cloud 
radar is developed to provide a more 
quantitative explanation of the observed 
water vapor evolution in the upper 
troposphere. In the selected case, regions 
with largest ice supersaturation ratios and 
deposition rates located were found near 
cloud top where plenty of small particles are 
present. This suggests a nucleation process 
in action. Moreover, a small particle regime 
is particularly sensitivie to specification of 
the deposition coefficient in terms of 
deposition rate and latent heat release. This 
implies that the cirrus dynamics and its 


evolution time scale may be closely linked 
to the deposition coefficient. Therefore, 
more studies to reduce the uncertainties in 
the deposition coefficient are warranted. 
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Figure 1. Flowchart. 
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Figure 3. Height-time sections of the derived 
D . , N t , and x • ln this calculation, v =2, 

J3 ; = 1, particles are columnar. 
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Figure 4. Heating rates due to latent heat 
release. 


Figure 2. Height-time sections of the RL 
extinction coefficient, radar reflectivity, and 
RHI. In the last panel, the black, white solid, 
and white dash-dotted curves indicate the 
cloud boundaries, the 30% error contour, and 
the 10% error contour, respectively. 



















